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= Royleigh  dissipation  function  for  structiaral  damping 


= modulus  of  rigidity. 


h = thickness  of  shell  wall 

Jq#  - Bessel  functions  of  real  argument 


[See  Equation  55] 


[see  Equation  56] 
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= natural  frequency 
= amplitude  of  propeller  force 

= generalized  force  due  to  fluid  reaction  on  end  caps 

[Equation  13 ] 

= generalized  force  due  to  fluid  reaction  on  cylindrical 
surface  [Equation  8] 

= generalized  propeller  force  [Equation  15 1 

= area  of  cylindrical  surface 
= cross  sectional  area 

I 

= thickness  of  end  cap 

= kinetic  energy  of  shell  cap  system 

» . longitudinal  displacement 

= longitudinal  generalized  displacement 

= potential  energy  of  deformation  of  shell-cap  system 

= radial  displacement  of  shell 

= radial  generalized  displacement 

= longitudinal  coordinate  ■ • 

= reactive  component  of  cap  acoustic  impedance  [Equation  34] 

= reactive  component  of  shell  acoustic  impedance  [Eqtiations 

55  - 57] 

= Bessel  Functions  of  real  argument 
acoustic  Impedance  of  end  capo 
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2/tm  for  m odd 

0 for  m even 


[See  Equation  43] 


[See  Eqiiation  6la] 


logarithmic  decrement  for  structural  damping 
Euler's  constant 

acoustic  Impedance  of  cylindrical  shell  surface  In 
water  [Equation  7] 

resistive  component  of  acoustic  cap  impedance  [Equation  34 

resistive  component  of  acoustic  shell  impedance 

[Equations  55-571 

angular  coordinate  of  shell 

mira 

& 

wave  length  of  sound  wave  in  water 
Poisson's  ratio 

mass  density  of  shell  material 
mass  density  of  water 
phase  angle  of  ra  the  mode 

^ ■ [See  Equation  55] 


ratr 

1^ 

& 

V P 

frequency  at  resonance  forcing  frequency  due  to  propeller 
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NOMENCLATURE 


a 


A' 

A. 


A,,  - 


mean  radius  of  shell 


X 

irphjSa  + 2p  c ira®  ~ 

^ ^0  0 CD 


2irphia  + 2irp  c mo  £ 
® ^ CD  2 


area  of  end  cap 


[Equation  42] 
[Equation  42] 


longitudinal  displacement  amplitude 


real  and  imaginary  parts  of  A 


B 


B 


2p  c^ira  0^  + K 
*^o  o c 


£,  — 

2irap  c 0 t:  + K 
*^o  o mo  2 


sound  velocity  in  water 


■ Vf 


■V 


E 

2p(l+v) 


= radial  displacemeni  amplitude 


C , C = 
r i 


real  and  imaginary  parts  of  C 


[Equation  43] 
[Equation  43] 


E 


Young's  modulus  of  shell  material 
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A SIMPLE  MATHEMATICAL  MODEL  FOR  AXIALLY  SYMMETRIC  MOTIONS 
OF  SUBMARINE  HULLS  DUE  TO  PROPELLER  EXCITATION 


I.  INTRODUCTION 


I 


A program  Is  cvirrently  being  pursuec^by  HYDRONAUTICS,  In- 
corporated^^ estimate  the  relative  degree  of  sound  radiation 
from  various  pulsating  sources  on. a submarine  hull.  This  particu- 
lar study  was  initiated  to  compute  the  approximate  effect  that 
the  pulsating  propeller  force  could  have  on  hulls  of  various 
diraenslonajjj 

Although  there  may  exist  more  accurate  mathematical  models 
for  describing  this  phenomenon  than  the  one  employed  here,  there 
do  not  exist,  to  the  writer's  knowledge,  any  closed  form  formulas. 
It  is  the  purpose  of  this  study  to  obtain  approximate  formulas 
which  can  be  used  to  assess  the  relative  importance  of  the  physical 
parameters  without  the  use  of  complicated  digital  programs. 
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II.  THEORY 

A.  Dlaplacemont  Components 

The  hull  will  be  approximated  by  an  elastic  cylindrical 
shell  with  rigid  circular  and  caps  (see  Figure  l).  The  rigid  end 
caps  may  move  longitudinally  as  a result  of  the  hull  shell  vibrat- 
ing in  a longitudinal  mode  of  vibration  but  the  caps  prevent  any 
radial  motion  of  the  shell  at  the  ends.  Only  the  axially  symmetric 
motions  of  the  system  will  be  considered  here  and  it  will  be  as- 
sumed that  when  the  shell  vibrates  in  a single  mode,  the  displace- 
ments can  be  approximated  as  follows: 

u = U(t)  cos  [1] 

. mTTX 

w = W(t ) sin  — 

In  equation  [Ij  m is  the  number  of  axial  half  waves  in  the  vibra- 
tion pattern.  For  the  first  mode  ra  = 1 and  it  will  be  seen  later 
that  the  hull  could  vibrate  with  displacements  that  are  primarily 
longitudinal  with  a small  radial  component  or  primarily  radial 
with  a small  longitudinal  component.  The  former  will  be  the  fun- 
damental longitudinal  mode  of  the  shell  and  will  be  of  relatively 
low  frequency  compared  to  the  latter  which  is  the  fundamental 
radial  mode. 

B.  Kinetic  Energy  and  Potential  Energy  of  Shell  Deformation 

Under  the  assumptions  of  axially  symmetric  motions  the  po- 
tential energy  of  deformation  for  the  m th  mode  of  the  shell  can 
be  vjrltten  as  follows 

V — [U^X^+W^+avTviTvJ  + [2] 

4a(l-v®)  12a=^ 


W.  Arnold  and  G.  B,  V/arburton,  Proc.  Roy.  Soc.  London,  A 197, 
23d-23ix  (194CQ. 
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and  the  kinetic  energy  is 

f = ZEgfa  + w*  ]*  . [31 

In  the  above  equations  7^  = . 

In  deriving  the  expression  for  the  potential  energy  the  assump- 
tions of  thin  shell  theory  have  been  employed.  For  a more  detailed 
account  of  these  assumptions  the  reader  is  referred  to  the  above 

reference^ , 

C.  Generalized  Coordinates  and  Generalised  Forces 

The  generalised  coordinates  of  the  system  will  be  chosen  as 
the  independent  displacements  U and  W.  The  generalised  forces 
corresponding  to  changes  in  these  displacements  will  be  associated 
with  the  following; 

1.  The  water  pressure  on  the  cylindrical  surface  arising 
from  radial  displacement  of  the  cylindrical  surface 

2.  The  water  pressure  on  the  end  caps  arising  from  longi- 
tudinal displacement  at  the  ends 

3.  The  force  on  the  hull  due  to  the  pulsating  propeller 

action. 


In  the  analysis  presented  here  the  following  assumptions 
will  be  made  vjlth  respect  to  the  water  pressures: 

1.  Tlie  pressure  on  the  ends  due  to  radial  motion  of  the 
cylindrical  surface  will  be  neglected  (this  assumption  should  be 
valid  for  primarily  longitudinal  motions) 

2.  The  pressure  on  one  end  due  to  motion  of  the  other  end 
will  bo  neglected  (this  assumption  should  be  valid  for  long  shells 


^The  end  cap  kinetic  energy  is  neglected  in  this  study.  This  is 
discussed  at  the  end  of  the  report. 
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where  the  ends  are  far  apart) 

3.  The  pressure  on  the  cylindrical  surface  due  to  motion 
of  the  ends  will  be  neglected. 

Let  the  generalized  force  associated  with  the  radial  pressure  (due 
to  radial  displacement  w)  be  then  the  virtual  work  due  to  a 

change  in  the  generalized  coordinate  ¥ is  6¥.  If  Is  the 

pressure  on  the  cylindrical  surface  due  to  radial  motion  then  the 
work  done  In  a displacement  6w  is 

/ p 6wds  (s  = area  of  cylindrical  surface) 
s 


8w  = 8¥  sin  ^ 


Therefore 


_ j£  j2rr  ^ m|2i  ad0dx 
0 0 ^ 


The  radial  pressure  p^  will  be  taken  equal  to  the  value  for  an  in- 
finitely long  cylinder  the  surface  of  which  has  the  radial  dis- 
placement proportional  to  sin  . The  value  for  this  pressure 
is  given  in  a previous  reference^  and  can  be  written  as  follows: 

p = -pc  W C (ka)  sin  — 7- 
^00  ^mo  ' m ' £ 

In  formula  [5]  is  the  density  of  the  medium,  c^  Is  the  sound 
velocity  In  the  medium  and  is  the  acoustic  impedance. 


'-M.  C.  Junger,  J.  Acouot.  See.  Am.,  25,  40-4?  (1953) 
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where  co  is  the  frequency  of  the  vibration.  The  acoustic  impedance 
can  be  written  in  terms  of  its  real  and  imaginary  parts  as  follows; 

C = 0 (k  a)  + i X (k  a)  [ 

^mo  mo  m ' mo'  m ' 

where  B is  the  resistive  component  (associated  with  radiation 
mo 

damping)  and  is  the  reactive  component  (associated  with  vir- 
tual mass). 

The  generalized  force  given  by  [4]  can  then  be  written  as 


follows ; 


o = > p c W C sln^  23^.  adOdx 


Q,.  = - 27rap  c W ~ ^ . [( 

^ ^o  o 2 mo 

Let  Qy  be  the  generalized  force  due  to  longitudinal  motion  of  one 
circular  end  cap,  then 

Q,.  6 U ^ - / p..  [6u]  dA  (A  = area  at  end  cap)  [i 

^ ^ x^o,£  ° ^ 

c 

In  equation  [9]  Py  Is  the  water  pressure  due  to  longitudinal  motion 
of  the  end  cap. 


Therefore 


0^.  ■-  - / p.4A  ==  P (force  on  end  cap), 
u . u c c 
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But 


= - Z , 
c c c 


[1< 


where  Is  the  acoustic  Impedance  of  the  circular  piston  end  cap 
moving  against  the  water. 

Thus 


Qy  = - Z-  U 


c c 

The  acoustic  impedance  of  the  end  cap  can  be  written  as  follovjs: 

Z = tra®  p c (9  + iX  ) 

c o ' c c' 

where  0 is  the  resistive  component  and  X is  the  reactive  compo- 
0 c 

nent.  The  values  for  0 and  X will  be  taken  as  the  values  for  a 

c c 

circular  piston  radiating  from  one  face  into  an  acoustic  medium. 
The  values  for  these  impedances  are  given  in  a recent  reference.® 
The  generalized  force  due  to  both  end  caps  v^ill  be 


[li 


Q, 


u 


20,^  r.  - 2p  c wa^  {e  + iX  ) U 
U ^0  0 c c' 


. -'1 

[131 


Lastly,  let  the  generalized  force  due  to  the  propeller  force  be 


Q , then 

Jr 


Thus 


= P = PBU 


= P = p e 
p o 


iO)t 


ui 


[15l 


D.  Equations  of  Motion 

Let  F be  the  dissipation  function  associated  with  structural 
damping,  then  the  Lagrange’s  equations  of/motlon  for  the  shell-end 
cap  system  will  be 


lest  Available  Copy 


, Hnnish,  IIRL  Report  U.  S.  Naval  Research  Laboratory, 

VJashington,  D.  C,,  Oct.  24,  i960. 


CONFIDENTIAL 


CONFIIffiNTlAL 


IKTOONAUTICS,  Incorporated 


-7- 


3t  = 


aw  aw  aw  ' ^ 


U61 


hij  Su 

where 

P = i K (U^  + W®)  [17: 

In  which  K is  proportional  to  the  structural  damping  force). 

Substituting  the  values  for  the  potential  energy,  the  kinetic 
energy,  the  dissipation  function,  and  the  generalised  forces,  the 
following  equations  are  obtained 

"irEhi  r 2 . TT  .*»  ^ j tr  \t»  . -r.  _iCOt 


^phj8aU  + 


2a(l-v=) 


[ + vxw]  + K il  = -2p  c 7ra®(©^  + iX  )U  + P e‘ 

Q \ 0 0 C X . o 


iirphiaW  + — + v>,u  + W]+  K W = -2Trap  c W 4 C 

" •^002  mo 


2a(l-v2) 


12a‘ 


[18] 


Assuming  that  the  motions  will  be  harmonic  in  time  the  equations 
can  be  rewritten  as  follows 


X 

c 


iirph;eaU+2p  cto®  U(— )+2p  Tra*0  U+K  U + 
^ ^00  'co  ^ ^o  c 


— — - [A*U+v>vW]=  P 

2a(l-v2) 


|TTphXaW+2irap^cJw(-~)  | + 2irap^c JlQ^^  ~ + — — ~-~['vJ+v?^u+  >v'*VJ] 


2 ^*0  o mo  2 2a(l-v2)' 


or  _ 

..  . 2p  c Tra=0  + K 
U+U  — 


irEhX 


+ K W = 0 


[7i^U  + v?.W3 


[19] 


P e 
o 


lent 


X X 

'irirphia+2p  c ira^  “ ^a^  x — - — 2 c 


o o 


0) 


2a(l-v®)  -|-Trphia+2p^c^ira^  ~ i7rphia+2p^c^Tra‘ 


c 

03  ■?, 


[STOP  0^9  I + icl 

W+Vj  — 


irEhj? 


uZ 

[W+v?^U+  Wl 

12a® 


X 


0 


^-irph£a+2irap^c^-^  I 2a{l-v2)  J-„.phia+2Trap^c^  | 


[20] 
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Before  going  further,  it  is  of  Interest  to  note  some  well  known 
special  cases  which  may  be  used  as  a basis  with  which  to  compare 
the  more  general  results  to  be  obtained  later.  The  first  case  is 
the  vacuum  uncoupled  longitudinal  vibrations  of  the  shell.  For 
this  case  the  equations  of  motion  reduce  to 

p 

P I = _°! [21] 

a(i-v^)  TTphia  ^phia 

The  natural  frequency  is  given  approximately  by 


P * ^ (neglecting  ) 

» V p X-v® 


[22] 


This  is  the  natural  frequency  of  a free-free  bar  of  length  i. 

The  second  special  case  is  the  in  vacuum  uncoupled  radial 
vibrations.  For  this  case  the  eqxiations  of  motion  reduce  to 


W + 

a(l-v2  j 


W+  — ?v*w| 

L 12a^  J 


Trphi^a 


0 


[23] 


The  natural  frequency  is  given  by 


]l  P(l-v=)  ^ V 


1 + 


12a‘ 


[24] 


As  Tv  -♦  0 this  reduces  to  the  well  known  formula  for  the  radial  vi- 
brations of  a shell.  Formula  [24]  Includes  a thiclcness  effect  term 

— — in  view  of  the  fact  that  bending  vjas  considered  in  the  ex- 

12a^ 

pression  for  the  potential  energy. 
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The  structural  damping  constant  K is  as  follows: 

K = ES£MS_  [25] 

where  p Is  the  natural  frequency 

6 is  the  logarithmic  decrement  for  structural  damping 
phia  are  as  given  before. 


E,  Solution  of  Equations 

Now  going  back  to  the  general  equations  and  letting 

„ -r-  lot  ,,  — iOJt 

U = Ae  , W = Ce 

the  following  equations  are  then  obtained  for  determining  the  am- 
plitudes A and  C 


A + icoA 


irEh/ 


1 s 

'2-7rph^a+2p^c^ira  ^ 


[ A+v>>  C ] 


2a(l-v®)  5TTph^a-r2p  c Tra®  — ^ 

^ ^ 0 Q CO 


X 

-2-7rph;2a+2p.  c ~ 

^ O O) 


[2Trap  c G :r  + K ] 
a . TT  ^o  o mo  2 

-0)  C+1cdC  r: + 


TTEh>e 


[ C+VXA  + — 


12a‘ 


C] 


[ 26a] 


^irphia+27rap  c 


mo  j2 


X 


= 0 [S5b] 


O O CD  2 


2a(l-v  ) ^ph;8a+27rap  c 


mo  £ 


O O CD  2 


Collecting  terms,  vie  obtain 


t!  TrEhjg 
A t 


j2a(l-v^)  ^ 

L 


ioj(2p  c Tra^G  +K) 

' o c ^ 2 

+ T; CD 


X . X 

•2-7Tph^a+2p  c Tra^  — ^ •-2irph;0a+2p  c — 

^ ^0  0 CD  '^0  0 CD 


TrEhi 


2a(l-v^) 


VA 


2T^ph^a-^2p  c ira^  — 
^ *^0  0 CD 


X 

■}:TrpbJa.+2p  c — 

' o 0 o 


[27a] 
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■rrEhjJ 


v?s 


,2a  (1-v*)  ^phia+2irap^c^-^  ^ 


2-^ 


+ C 


irEhi 


■(1+ 


h* 


2a(l-v^)  12a  ‘ 


^ ^ io)(2Trap  c 6 4 + K 

' ^o  o rao  2 

1 -ir  rn  -I  ""  ' — I I r - ■ T 


S. 

5irhipa+2Trap^c^  “ iTrphia+2Trap^c^' 


X 


-05' 


= 0 


[27b] 


mo  ^ 

CO  2 


Equation  [27]  constitute  a set  6f  two  complex  simultaneous  alge- 
braic equations  which  are  to  be  solved  for  the  complex  amplitudes 
A and  Co  The  solutions  will  be  of  the  following  form 


A = A + i 

r i 

C = C^  + 1 C^ 

where  the  r subscript  refers  to  the  real  part  and  the  1 subscript 
refers  to  the  imaginary  part.  The  displacements  themselves  will 
. then  be  equal  to 


[28] 


wh(?re 


u - A^®  + A^^  cos  cos 

w = sin  cos 

-1  ^i 

♦ = tan  

m ft 


(cot  + b ) 
' m*^ 

(cot  + 4>  ) 


[29] 


r 

is  the  phase  angle  between  the  driving  force  and  the  displace- 
ments and  arises  because  of  the  radiation  damping  and  structural  damp- 
ing in  the  system. 
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P.  Longitudinal  Motions 

The  motions  that  are  of  prime  importance  as  a result  of  pro- 
peller excitation  are  the  primarily  longitudinal  ones.  These  un- 
coupled motions  in  the  water  will  now  be  studied  in  more  detail. 

The  equation  of  motion  for  uncoupled  longitudinal  motions  in  the 
water  is 


..  U [2p  c ira*©^  + K] 
tl  + LS-O £ ^ + 


TrEhX 


U 


■ X 

c I,  /,  ..2  V 1 . . _ a c 


iirphia+2p^c^Tra2  ~ 4a{l-v^)  ^ph£a+2p^c^7ra 

Icot 


CO 


P e" 
o 


27rphia+2p  c 

^ *^0  0 CO 


[303 


Assuming  U = A e^^^,  we  obtain 


^ ico(2p^c  vtL^0^  + K)  _ . 

Q 0 ' ^ jrgh^ 


2?v‘ 


^phia+2p^c  ira®  -2. 

^ ^0  0 CO 


>8  C 


4a(l-v  ) ^7rphi-a+2p  c ira* 

. ' ^ *^0  0 CO 


X 

•§irph-2a+2p  c Tra^ 

*^0  0 CO 


[31] 


P 

o 


•g-TTpiiia  + 2p  c Tra®  ~ 

_V *^0  0 CO 


[32] 

ico(2p^c  7ra®0^  + K) 

■ 

4a(l-v®)  ■^•7rph£a+2p^c^Tra®  ^ 

- 

••  U) 

X 

■i-rrph-fia+2p^c  im® 

*^0  0 , CO 

Tl^e  natural  frequency  of  these  longitudinal  motions  in  the  water 
lo  therefore  given  by  the  following  frequency  equation 
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irEhi 


2'K 


a 


- CD*- 


0 


4a(l-v^)  iirphj8a+2p^c^Tra®  ~ 


l33] 


The  values  of  o>  which  satisfy  [33]  are  the  natural  frequencies  of 
the  system.  In  general  co  cannot  be  obtained  directly  since  is 

a function  of  co.  For  relatively  low  frequencies,  however,  the  re- 
sistive and  reactive  impedance  of  an  unbaffled  circular  piston 
radiating  from  one  side  into  an  Infinite  medium  can  be  approximated 
as  follows®: 


7 i ka  [34] 

where  k = ~ . 

X c 

o o 

Under  these  circumstances  the  frequency  equation  becomes 
vEhi  2X® 

; — - CD®  = 0 

4a(l-v2)  |Trphjga+7ra3p^ 

The  natural  frequency  can  then  be  obtained  in  closed  form  and  is 
as  follows; 


[35] 

3est  AvailBblB  Copy 


irEh^  2^2 


~i-TTphjCa+ira3p 


or  neglecting 


CO 


1+  2 


a®  P< 

hir 


[36] 
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So  the  natural  frequency  of  the  low  frequency  longitudinal  notions 
in  water  is  equal  to  the  .in  vacu\iia  frequency  reduced  by  a factor 
which  depends  on  the  radius  of  the  shell,  the  thickness  of  the 
shell  wall  and  the  shell  length.  It  is  interesting  to  note  that 
the  factor  in  the  denominator  of  equation  [36]  can  be  written  as 


aa^p, 


hip  27rahip 


3M 


M 


sphere 


[37] 


shell 


where  Is  the  mass  of  a sphere  of  water  which  has  the  sane 

radius  as  the  shell  and  M . .is  the  mass  of  the  shell  without  the 

Gxl0XX 

end  caps.  When  the  shell  is  oscillating  at  its  natural  frequency 
the  amplitude  of  oscillation  is  found  from  equation  [32]  and  is  as 
follows; 


m (2p  c ra®0^+K) 
n *^0  0 c 


[38] 


Using  the  low  frequency  approximation  for  0 the  following  expres- 

0 

Sion  is  obtained I 

P 

A|  = -2 


2 2 

/«  2 CD  a 

cD(2p  c tra  + 

' ^o  o 


coSphia 

2 


or 


in  = 


p Tra  CD 

icD^phia(-2.  — p ^ +6) 

n*^  ' p c hi  ' 

• r\ 


[39a] 


The  propeller  force  is  about  10  percent  of  the  xneai^  thr-uct^  the 
mean  thrust  is  equal  to  the  drag  and  the  drag  is  civin  approximately 
by  the  following  formulas 
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D = i 

where 

D = mean  drag 

Cjj=  drag  coefficient  (about  .05  for  modern  submarines) 


V = velocity  (15-30  knots) 
S = cross  sectional  area 


G.  Coupled  Longltudlnal.Radlal  Motions 

Going  back  to  the  more  general  equations  [2?]  the  second  of 
the  equations  gives  the  following  relation  between  the  radial  dis- 
placement amplitude  C and  the  longitudinal  amplitude  A for  a given 
longitudinal  driving  force 


J 


rrEhje 


2VA 


C = - 


o O CD  2 


_1x 

no  £ J 


TEhi 


4a(l-v°) 


(2  + 2 


— 'X'*)  + Ico  (2irap  c 9 4 + K) 

a ' ^o  o mo  2 ' 


12a 


iirphia  + Pimp  c 4 

^ *^0  o CD  2 


- 0) 


640] 


Therefore  substitution  into  the  first  of  equations  [27]  gives 

P 


A = 


a2X 


B 


A ^ A - - 


o2vx 

a2v>v 

La  J 

LalJ 

where 


a (2+2  X^) 


12a‘ 


X 

A = ir7Tphj5a+2p  c Tra^ 

^ ^0  0 CD 


^ . B'  a 

+ ICD  ~ - CD 


X . 

A ' = ■girph^a+Slirap  c — ^ -r- 
'^O  O CD  2 


[41] 


[42] 
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^ — irMiL.  B = 2p  0 ira^e  + K B'  = 2irap  o 4 +k 

4a(l-v=)  »'o  o 0 *-0  o mo  2 

K Is  the  structural  damping  parameter  as  given  before. 
Simplifying 


[a2>v®+iajB-Aai®]  - 


[a2VA]  [a2v?v] 

[a (2+2  — 7^^)  + ioiB'  - A'oi®] 


[44]j 


Now  divide  the  denominator  of  the  above  equation  into  its  real 
and  imaginary  parts  as  follows! 


Denominator  « 


[a2^®+icDB-Ao)®]  [a(2+2  )+lo3B' -A’cd®] -[a2V>]® 

12a® 


[a(2+2  — 7v^)  + icoB'  - A'co®] 


K = a2X®  - Acs® 


1, 9 

K'  = a(2+2  ~ — ?v^)  - A 'CO® 
12a® 


Then 


[451 


m- 


[46]? 


Denominator  = 


;k+1coB]  [K'+icoB']  - [a2v?v; 
K’  + icoB' 


[47] 


= (K  - 


[a2v?v]®K' 


[K'®+co®B'® 


+lco  ( B+ 


[g2v?c]®B' 
K'  ®+co®b'  ® 


For  uncoupled  longitudinal  motions  set  v o and  the  equation 
reduces  to  the  uncoupled  longitudinal  case  obtained  previously. 


[48] 
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The  natural  frequency  for  coupled  motions  is  given  by  the  solution 
to  the  following  frequency  equation: 


K - K' 

[K'®+05®  B'®] 


[491 


and  the  resonant  amplitude  of  the  longitudinal  component 
is  given  by 


P 

o 


CD  ( B+ 


ta2v^]^B'  1 


of  motl'on 


[50] 


where  cd  Is  the  natural  frequency.  For  small  values  of  cd  the  term 
cd®B’®  will  be  considerably  smaller  than  K'®.  Under  the  assumption 
that  o)®B’®  « K'®  the  frequency  equation  becomes 

KK'  = [a2v>]®  [51] 

At  low  frequencies  it  can  be  shown  (and  will  be  shown  later)  that 
A and  A'  are  independent  of  cd.  Therefore  the  frequency  equation 
reduces  to 

la29^®][o(2+2>^  7s^)]-cD®([A][a(2+2-^X^)]+[A'][a2X®]  ) +AA'co^ 

12a®  I 12a®  i 


x=  [a2v>)®  [52] 

If  we  focus  our  attention  on  modes  of  relatively  long  wave  length 
where  X is  small  and  if  we  neglect  v®  compared  to  unity  in  the 
above  equation  then  the  shell  end  cap  system  has  the. .following  two 
sets  of  natural  frequencies: 


CD* 


2aV 


[53] 


These  frequencies  arc  exactly  the  uncoupled  radial  and  longitudinal 
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frequencies  in  water  respectively.  The  lower  one,o>  , is  the  long!- 
tudlnal  and  Is  the  one  of  greatest  practical  Importance  in  this 
study. 

In  order  to  simplify  the  theory  further  some  rather  far 
reaching  assumptions  have  to  be  made  concerning  the  values  of  B, 

B*,  and  K' . Therefore  the  following  section  will  be  devoted  to 
the  simplification  of  these  expressions  and  the  Justification  for 
this. 


If  it  is  assTuned  that  the  impedance  offered  by  radial  motions 
can  be  approximated  by  the  Impe^ancec  cifr  an  infinitely  long  cylinder 
in  an  infinite  acoustic  medium  whose  vibration  pattern  is  identical 
to  the  finite  shell  (see  recent  reference^)  then 

e I 

B'  = 2irap^c^  [54.3 


If  Y c 


K'  * a (2+2 
r/°o<" 


12a‘ 


A ' = ^phj2a+2Trap  c 


mo  I 


O O CD 


e =0  j 
mo 


X = 
mo 


a) 

(k;a)Kjk;a) 


where 


o 


P = 


2 p 

1+2 

h£  p 


Greenspon, 
August  i960. 


Jour.  Acoust.  Soc.  Am.,  Vol.  32,  No.  8,1017-1025, 


CONFIDENTIAL 


CONFIDENTIAL 


HYDRONAUTICS, ^Incorporated 


-18- 


niT 


for  resonant  longitudinal  vibrations  p = o>  P 


■/r 


k'a  = 
m 


M c /c^[J^{k  a)j  (k  a)+Y„(k  a)Y  (k  a)3 
Tr  o o m 1 m o m 1 m 


mo 


(K^)  I t-J  (k  a)l®+[Y  (k  a)]^ 


m 


I 


1 ' n 


1 ' m 


[55] 


0 


2 A Y c /c 
r'  o 


mo 


whore  )P^  - 1 = XK« 


K* 


p®  - 1 


For  small  k^^a  (small  ?y)  and  Y c > 1 

>^Y(c  /c  )Y  (k  a) 

„ ' ir  o'  o'  m ' 


mo 


mo 


Tk  a)Y  (k  a) 

gM(c/c^) 

ir(h„a)^  V {k  a) 
m 1 m 


[56J 


where 


k a 
m 


Y (k  a)  ~ — (j£n  — ^ + t)  7 = Euler's  Constant 
o in  ^ 


TT 


Y (k  a)  = ~ — 
X m vk_a 


ra 


r^' 
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So 


(te  V* 

®mo 


HiflU.UU  J(IU,II.  i|..  . 


X * 

(therefore  a (&n  + 7)  so  that  A’  Is  Independent  of  the 

mode  and  frequency  for  small  (a). 

Thus  (neglecting  structural  damping) 


B*  = 


2TOp^Co  I > (c/c^)P  i 


K'  = a (2+2  )-<u^ 

12a" 


lirphia  + 


2ira®p^c^>v  ^°e/®o  I 

____ 


Eh£ 

a 


2ir 


A(l-v2) 


- I 27T  ^ (£n  ^ + 7) 


For  email  X the  last  two  terras  r\fill  he  neglected^'  leaving 

j^i  ~ Ehf  ^ 27t  j 

^ I'.(l-v2) 

Thus  inserting  the  above  relations  in  Equation  [50I 


A 


03  < 2p  c Tra*"©  -I 
■ o c 


(a2vx)227Tap^c^  | X (c^c^  ~ 


"Eh<^ ' 


27T 

4(l-v2) 


Duo  to  the  behavior  of  the  logarithm  this  assumption  may  be  in 
considerable  error'’  for  some  geometries. 
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Slmpliryinc  and  su'ootituting  the  value  for  0^ 


o)/p^c  Tragi'  +v^’K^2Trap 

o 


r A1 

o^o2  c ' 2J 
o 


[60] 


§pca^~p(~P  + v^m  ttA 
2 *^o  o c 1 c j 

r\  r\  * 


[61] 


Instead  of  employing  the  resistive  impedance  of  an  infinite  cylin- 
der let  us  use  Robey' s®  results  for  the  resistive  impedance  of  a 
uniformly  pulsating  ring  and  average  the  displacement  value  for  the 

given  mode.  'The  B'  can  then  be  written  (for  small  k ) 

a 


[6la] 


where  £ is  the  length  of  the  cylinder  and  k = 


r2r  . rmrx  . , * 

J J sin  -7—  ox  a d 0 

0 0 

/ / adOdx 


2Tra  — (-cos 

■n^>rr  ' 


m7T 

27raf 


S, (0 


— 2 

,3  = — if  m is  odd 
rav 


3=0  if  rc  is  even. 


Cubstituting  this  into  the  equation  for  A,  we  obtain 


1 ^a-xa  H-a 


o o 2 


} ( l+(3  2 ^ 


[62] 
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or 


[63] 


CD  ( f a®7^*  ~ F*  ) (l+2v*  p ra^Ti^ 

c^ 


o'^o  2 » 


We  will  use  this  expression  instead  of  [61J.  The  amplitude  of 
radial  motion  at  longitudinal  resonance  can  be  written 


r q2v>>  , 

»*  A I J 


C = - A 


or 


(2+2 


a 


12a‘ 


. 4 B'  s 

+ iCD  7-; CD 

A ' 


[64] 


C = - 


P [a2vx] 
o 


[K+1cdB][K'+1cdB' ] - [2av>,]  = 


[65] 


neglecting  co  BB'  compared  to  the  other  terms  we  obtain 
C 


P^  [a2V7v] 
o 


[KK'-{2VaA)“]  + ico(K'BfB'K) 

At  resonance  KK'  = (2Va?\)*  (see  equation  [51])  therefore 
radial  amplitude  at  longitudinal  i^sonance  is 


C » - 


P^[a2v>v] 

o 


cd(k'B  + B'K) 

Substituting  the  value  of  the  longitudinal  resonance  frequency, 
following  relation  is  obtained  for  C 


C = - 


P^[  2v?^a] 
o 


u)(K'  2ira^p  c^ 
*^00 


CO 


4c= 

o 


_ c® 

a^+K  J3  p c tt£-K^  — P®  i ) 
'^'^00  2 ' 


[66] 

the 

[67] 

the 

[68] 


lest 
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C = - 


[2VaX] 


[69] 


r I a 


The  volume  change  of  the  coupled  motions  at  longitudinal  resonance 
can  bo  written 


AV{t)  = wdxadO  + (U  U .) 

Q Q 


[70] 


— 2TraC(-cos  rar+cos  0)+  (U 
mu  ' x=£  x=0'^  cap>e 


i(o3t+^)  [71] 


where 


♦ = phase  angle. 

Therefore  only  for  odd  m do  we  obtain  a net  volxinc  change  due  to 
radial  motion.  For  modes  with  m even  there  is  no  volume  change  at 
all  because  the  ends  are  moving  in  phase  with  the  same  displace- 
ment. The  volume  change  can  be  written  as  follows: 

For  m odd 


AV 


P^[2v>a] 


CU(P  c 
yo  o 


I 


P 2Tra‘ 
o 


„2,a  e ,.3 


^{Po°o  2 "T  /3 


o o 2 


o 


3o‘ 


s\ 
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For  m even 

AV  --  0 [75] 
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discussion 

Approximate  closed  formula  relations  have  been  derived  in 
this  report  for  the  eouplod  longitudinal  - radial  motions  of  a 
cylindrical  shell  - end  cap  configuration  which  is  intended  to 
model  in  a rough  manner  a submarine  hull  undergoing  axially  sym- 
metric motions,  A niimber  of  rather  drastic  simplifying  assumptions 
have  been  made  in  order  to  obtain  the  closed  form  relations.  The 
rewards  of  such  a simplified  analysis  are  that  the  relative  effects 
of  different  jsarameters  can  be  assessed  and  problem  areas  can  be 
defined  without  resorting,  to  detailed  computer  analyses. 

This  simplified  approach  can  be  used  to  point  out  salient 
features  but  it  is  believed  that  the  assumptions  bear  some  further 
investigation.  For  the  future  it  is  believed  that  the  following 
items  should  be  examined: 

1.  Investigate  the  accuracy  of  the  simplifications  by  carry- 
ing out  computer  calculations  with  the  unslrapllfied  equations  [26] 
and  [27] . 

2.  Compare  the  results  with  available  spheroidal  shell 
analyses , 

3.  Compare  the  relative  magnitude  and  frequency  of  sound 
radiated  by  the  hull  vjith  sound  radiated  by  other  sources  now  under 
investigation. 

4.  Consider  non-resonant  response  in  combinations  of  modes. 

5.  Compare  the  analysis  viith  available  experiments. 

6.  Determine  the  sound  field  in  more  detail. 

7.  Consider  the  effect  of  transverse  rings  and  bulldieads 
in  tiiG  shell  analysis. 

<j.  Consider  tfie  effects  of  structural  damping  in  more  detail. 

9.  T4a]ce  parametric  studies  using  existing  hulls  as  examples. 

10.  Jlstliaate  the  o of  the  various  modes. 

11.  listlmatc  the  effect  of  mutual  Impedance  between  the  ends 
a.'V.i  cylindrical  r^ur’faoc. 
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In  thlc  c:tudy  wo  have  found  that  volume  changes  are  associated 
vjlth  the  odd  m modes  only,  the  volume  change  being  given  by  [7^-]« 

The  even  m modes,  on  the  other  hand,  arc  not  associated  vjith  volume 
changes  and  essentially  act  as  a dipole  radiator  vjith  a source  at 
one  end  of  the  shell  and  a sink  at  the  other. 

If  Gtiffonors  and  rigid  bulldieads  were  added  this  should 
tend  to  decrease  the  radial  motion  and  allow  for  more  of  a net  vol- 
ume change  due  to  combined  longitudinal  motion  and  radial  motion. 

Supposing  bullcheads  were  inserted  which  were  allowed  to  ro- 
tate but  would  prevent  radial  motion  at  their  attachment  points. 

Tliis  would  probably  give  a radial  displacement  distribution  resembl- 
ing a continuous  sine  curve  with  negative  and  positive  crests.  If 
the  number  of  hills  and  valleys  of  radial  displacement  were  equal 
then  the  volxirae  change  to  due  radial  motion  would  cancel  leaving 


4 P 


AV  = 


m odd  only 


vjhere  cd  Is  the  natioral  frequency  of  longitudinal  hull  vibrations 
and  can  be  written  as  follows: 
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In  this  analysis  we  neglected  the  kinetic  energy  of  the  end 
caps.  This  tartn  would  be  i M U*  where  M =»  ira*pt  (mass  of  end  cap, 
t ■ thickness  of  cap).  The  ratio  of  longitudinal  kinetic  energy  of 
the  shell  to  the  kinetic  energy  of  the  two  end  caps  would  be 

h i 

W a 

We  limited  the  analysis  to  cases  where  / » a (h  will  be  of  the -same 
order  as  t). 

In  order  to  consider  end.  cap  kinetic  energy  the  total  kinetic 
energy  would  be 

T = + ) 

Therefore  we  merely  carry  the  analysis  through  with  the  longitudinal 
inertia  terra  increased  by 
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